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Colloidal 2D monolayers sliding in an optical lattice are of recent importance as a frictional sys¬ 
tem. In the general case when the monolayer and optical lattices are incommensurate, we predict 
two important novelties, one in the static equilibrium structure, the other in the frictional behavior 
under sliding. Structurally, realistic simulations show that the colloid layer should possess in full 
equilibrium a small misalignment rotation angle relative to the optical lattice, an effect so far unno¬ 
ticed but visible in some published experimental moire patterns. Under sliding, this misalignment 
has the effect of boosting the colloid monolayer friction by a considerable factor over the hypothet¬ 
ical aligned case discussed so far. A frictional increase of similar origin must generally affect other 
incommensurate adsorbed monolayers and contacts, to be sought out case by case. 


PACS numbers: 68.35.As,68.35.Gy,83.10.Rs,82.70.Dd 


The mutual sliding of crystalline lattices offers, despite 
its apparently academic nature, one of the basic plat¬ 
forms to understand the nanoscale and mesoscale fric¬ 
tional and adhesion phenomena [1]. In one of the fresh¬ 
est developments, Bohlein and collaborators [2] showed 
that the sliding of a 2D crystalline monolayer of colloidal 
particles in an optical lattice provides unexpected infor¬ 
mation on elementary tribological processes in crystalline 
sliding systems with ideally controlled commensurabili- 
ties. Given the scarcity of reliable and controllable fric¬ 
tional systems, it is hard to overestimate the importance 
of such model systems with full external control over all 
parameters including periodicity, coupling strengths, and 
applied forces. For this reason 2D monolayers in periodic 
lattices require a close theoretical study. In this letter 
we describe two main surprises, one structural and one 
frictional, which emerge from realistic molecular dynam¬ 
ics simulations. We show first of all that incommensu¬ 
rate colloid islands naturally develop in full equilibrium a 
small misalignment angle relative to the substrate. Sec¬ 
ond, sliding simulations demonstrate that the misaligned 
angular orientation increases significantly the dynamic 
friction with respect to the (hypothetical and unstable) 
aligned case. While both are important for colloidal 
monolayers, their potential impact extends in principle 
beyond the specific case, to a wider variety of systems 
where mutually incommensurate 2D lattices are brought 
in static and then in sliding contact. 

In colloid monolayers, the 2D density may vary from 
“underdense” (p < 1, where p = a//a c , with ai the spac¬ 
ing of the optical lattice, a c that of the colloid lattice, 
both triangular) to perfectly commensurate (p = 1, one 
particle per potential well), to “overdense” (p > 1), each 
with its specific sliding behavior. Both experiments [2] 
and theory [3, 4] indicated that commensurate (p = 1) 



FIG. 1. (Color online) (a) Schematic of a 2D colloidal lattice 
(black dots) aligned with an incommensurate triangular pe¬ 
riodic potential mimicking the optical lattice (blurred spots 
represent potential minima), (b) A misaligned configuration 
rotated by an angle 6. 


friction is large, dropping to much lower values for p ^ 1, 
where optical and colloid lattice are incommensurate. 
This drop reflects the great mobility of the pre-existing 
misfit dislocations, also called kinks or “solitons”. In 
these studies the two lattices, colloid and optical, were 
silently assumed to be geometrically aligned, prior to and 
during sliding. That assumption however is dangerous. 
A long-known theoretical result suggests [5], for example, 
that a harmonic monolayer subject to an incommensu¬ 
rate periodic potential of weak amplitude Uq may partly 
convert the misfit compressional stress to shear stress by 
an equilibrium geometric misalignment of the monolayer 
(see Fig. 1) through a small rotation angle 

f 1 + p 2 ( 1 + 25) \ ( ^ 

^ — (1) 

whose energy-lowering effect originates from a better in- 
terdigitation of the two lattices. Independent of [To, the 
rotation angle is nonzero in this approximation provided 
the transverse 2D sound velocity cp is sufficiently smaller 
than the longitudinal cp, precisely if 5 = (cl/ct) 2 — 1 > 
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p _1 . While this kind of rotated epitaxy has been ad¬ 
dressed experimentally [6, 7] and theoretically [8] for ad¬ 
sorbed rare-gas monolayers, its possible presence in col¬ 
loidal monolayers was so far unsuspected. More impor¬ 
tantly in the context of sliding friction, the tribological 
impact of an equilibrium geometrical misalignment is un¬ 
explored in any incommensurate system. The externally 
forced rotation turning a commensurate layer into incom¬ 
mensurate is well known to reduce friction, as exemplified 
by the sliding of graphene flakes on graphite [10, 11]. Dif¬ 
ferent as these two cases are, a possible expectation might 
be that the equilibrium geometry, alignment of commen¬ 
surate layers, or misalignment of incommensurate layers, 
should always exhibit a higher friction relative to forcedly 
rotated ones, since the optimal T = 0 geometry must in 
every case corresponds to a closer interdigitation of the 
two lattices. 

First, let us consider structural alignement. Using the 
same methods as in Ref. 3, we model the colloidal system 
in an optical lattice as a 2D monolayer of N p point-like 
classical particles, mutually repelling through a screened 
Coulomb potential U pp (r^-) while immersed in a static 2D 
triangular lattice potential Wfa) = Uow{ri) where w{ r) 
is a dimensionless periodic function of spacing ai ( as 
specified in Supplemental Material ( SM) ), and Uo is the 
amplitude (“corrugation”) parameter. The Hamiltonian 
is thus 


N v 




i= 1 


U 0 w(ri ) + I 


3 7 ^ 
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The particles are confined to the (x,y) plane, and sub¬ 
jected to either planar periodic boundary conditions 
(PBC) in a given area A, or alternatively to an additional 
Gaussian confining potential G(r ) = —exp(—r 2 /cr^), 
in which case they form an island with open boundary 
conditions (OBC). Temperature is generally set to zero, 
because the results are clearer and require less statistics 
in this limit. Finite temperature results are otherwise 
not essentially different, as shown along with details of 
optimization protocols in Supplemental Material {SM). 
Within the confinement region the 2D particles crys¬ 
tallize in a triangular lattice of mean spacing a c , with 
modulations induced by the periodic potential (in the 
OBC case there is also a confinement-induced variation 
between a dense center and a sparser periphery). Here 
and in the rest we shall focus for specificity on the un- 
derdense incommensurate case p = 3/(1 + \/5) — 0.927; 
the physics with different values of p ^ 1 including p > 1 
is qualitatively similar, as briefly discussed in SM. 

We determined, by means of careful energy minimiza¬ 
tion, the optimal T — 0 2D geometry of all N p particles 
for increasing corrugation strength; shown here are the 
results for Uo = 0.1 — 0.6, where the effects are particu¬ 
larly clear. The final, optimal geometry of the 2D colloid 
lattice initially aligned at 6 = 0° with the optical lat¬ 



FIG. 2. (Color online) Relative static energy of structure- 
optimized colloid islands (OBC) as a function of the trial rota¬ 
tion angle 0. (a) Total energy per particle E to t- (b) Periodic- 
potential (corrugation) contribution W to E to t . (c) Interpar¬ 
ticle interaction contribution U pp to E tot . Curves correspond 
to increasing corrugation amplitude Uo = 0.18 — 0.54. Ener¬ 
gies are measured relative to that of the colloidal monolayer 
at rest and at Uo = 0. Dashed line: ideal NM angle, Eq. (1), 
Onm ~ 2 . 6 °. 


tice axes is found to be misaligned, with a small rotation 
angle 0 opt — 2.3° in PBC calculations. This rotation real¬ 
izes a better interdigitation with the optical lattice, and 
occurs spontaneously during the simulation at the cost 
of creating the dislocations required by the PBC con¬ 
straints (see SM). More detailed energy minimizations 
were done in OBC, which do not have the same prob¬ 
lem. Since the 30,000 particle islands are too large to 
spontaneously rotate, we carried out simulations start¬ 
ing from a prearranged sequence of misalignment angles. 
The total energy minimum versus 6 confirms, as shown 
by Fig. 2, the structural misalignment angle at equilib¬ 
rium, with a magnitude in the range 3° < 0 opt < 5° 
(depending on the optical lattice strength Uo) generally 
somewhat larger than in PBC, where angular constraints 
hinder the misalignment. 
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FIG. 3. (Color online) Equilibrium configurations obtained 
for p ~ 0.83 and Uo = 6.3. (a) Unrotated 0 — 0°; (b) opti¬ 
mally rotated 6 opt = 7°. Dark/light colloids enjoy best/worse 
W. Only the central part of the island, optimized in OBC, is 
shown, (c) Experimental geometry for the same p, adapted 
from Ref. 12, where both the moire angles and spacing com¬ 
pare directly with (b) rather than (a). 


These structural results compare instructively with 
those expected from the harmonic model [5]. In our case 
a 2D phonon calculation for the monolayer yields a sound 
velocity ratio cl/ct = 1.806, larger than the theoretical 
threshold value (1 + p -1 ) 1 / 2 1.442. The correspond¬ 

ing theoretical misalignment #nm — 2.6° is in qualita¬ 
tive agreement with the more realistic simulation result. 
Figure 2b,c shows how the two pieces which compose 
the total energy, namely the periodic lattice energy part 
W = (W(ri)) controlled by the corrugation amplitude 
Uo and the interparticle interaction C/ pp = (77 pp (Gj))> be¬ 
have. Misalignment raises the interparticle energy, but 
that cost is overcompensated, by a factor 2, by a corru¬ 
gation energy gain. The incommensurate colloid static 
equilibrium structure is misaligned because that permits 
a better interdigitation with the optical lattice. 

Even if our predicted misalignments are small, their 
experimental existence is easily revealed, because the ro¬ 
tation angle 0 between two lattices is highly amplified by 
the moire pattern, which rotates relative to the periodic 
potential lattice by an angle ^ satisfying the geometric re- 
lation cos 6 = p -1 sin 2 ip + cos 'ip [l — p -2 sin 2 !p] [16]. 

As the moire superlattice rotates by ip, its spacing L also 
decreases [17] from its aligned value of about L=a c p/{ 1- 
p) to its rotated value of L'=a c /y/l + p -2 — 2p -1 cos 0. 
As an example, Fig. 3 reports the structures of the artifi¬ 
cially unrotated and of the optimally rotated ( 0 opt = 7°) 
configurations calculated for p ~ 0.83 and Uo = 6.3 (pa¬ 
rameters believed to be appropriate to experiments in 
Ref. 12) in comparison with one another and with the cor¬ 
responding experimental structural moire. Both the ori¬ 
entation and spacing of Fig. 3c agree with the 0 opt = 7° 
but not with the 0 opt = 0° pattern, proving that the 
misalignment was actually present in that experiment. 

The particle static displacements associated with the 
optical lattice potential are also enlightening. Figure 4 
shows the moire pattern of a small portion of the mono¬ 


FIG. 4. (Color online) Moire patterns of the particle mono¬ 
layer’s central region as obtained in OBC for p = 0.927, 
Uo = 0.27, (a) for 6 — 0° (moire angle 'ip = 0°) and (b) 
for the optimal 0 = 4° (moire angle ip ~ 15°). Each dot 
indicates a particle in the unrelaxed configuration, colored 
according to the local corrugation potential W{ r): dark for 
potential minima, bright for maxima. White arrows show the 
displacements of each particle from the ideal triangular lat¬ 
tice to the fully relaxed configuration, magnified 15 times. 
The compression-dilations at 6 — 0° are turned into largely 
shear, vortex-like displacements at 0 opt = 4°. 


layer island (p = 0.927, Uo = 0.27) for 6 = 0° and for 
#o P t = 4°, corresponding to a moire angle ip « 15°. 
Particle displacements, designated by arrows, change 
from radial longitudinal compression-dilations to mixed 
shear-longitudinal, vortex-like displacements upon opti¬ 
mal misalignment. A large 2D bulk modulus and a weak 
shear rigidity of the crystalline monolayer are crucial fac¬ 
tors increasing the extent of the shear distortions, there¬ 
fore enhancing the lattice misalignment. 

We come to our second point, i.e., the forced sliding 
of the particle monolayer over the periodic corrugation 
and the associated frictional losses. The shear distor¬ 
tions and the corresponding increased interdigitation at 
the optimal misalignment angle 0 opt are expected to af¬ 
fect the sliding of the particle lattice over the periodic 
potential. Sliding is realized by a flow of the soliton su¬ 
perstructure, accompanied by dissipation as part of the 
work goes into soliton-related time-dependent distortions 
of the 2D lattice. That work will change once the nature 
(longitudinal to shear), orientation (0° to -0), periodicity 
(L to L') change with 0 (0° to 0 opt ). We determine the 
magnitude of the expected friction change by simulating 
the overdamped sliding dynamics of the OBC island over 
a range of 0 values, so as to assess the frictional effect 
of misalignment near its optimal value. We applied an 
external driving force F^ acting on each particle, slowly 
varying to-and-fro as a function of time, mimicking the 
experimental drag force rjVd induced by a fluid of viscos¬ 
ity rj and slowly back and forth time-dependent speed 
Vd [2] (details in £714). Despite a nonzero torque, gen¬ 
erally present for all preset angles that differ from # opt 
(and from zero) the misalignment angle did not have the 
time to change appreciably in the course of the simula¬ 
tion. Under sliding, the frictional power dissipated per 
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FIG. 5. (Color online) Dissipated friction power pf r i c as a 
function of the trial misalignment angle 6 (Uo = 0.27). Three 
curves are reported corresponding to increasing values of the 
driving force Fd — 0.18, 0.27, 0.36. The inset shows pf r i C as 
a function of Fa for two values of 0. 


particle was calculated as [3] 

Pfric = F d -(v cm )-7?|(v cm )| 2 = (r?/AT p )E i (|v i -v cm | 2 ) (3) 

where v cm is the center-of-mass velocity, v* is the velocity 
of particle i, and brackets denote steady-state averages. 
Due to the confining envelope potential, the lattice spac¬ 
ing a c , close to constant in the central part, increases to¬ 
ward the periphery, where colloids also tend to be pinned 
by the corrugation. To address properties of mobile col¬ 
loids at a well-defined density, trajectories were analyzed, 
as also done in experiments, considering only particles 
belonging to the central portion of the island, in our case 
a square of size 80 x 80 a^. The duration of the sliding 
simulations was fixed by requiring a total center of mass 
displacement not smaller than Ax cm « 2 — 3a c . 

Figure 5, our main dynamical result, shows that fric¬ 
tion is increased by a very substantial factor by mis¬ 
alignment relative to alignement, reaching a maximum of 
about two at the optimal angle 0 opt , subsequently drop¬ 
ping for larger angles where the energy gain and static 
distortion magnitude also drop. The physical reason for 
the frictional peak at 0 opt can be further appreciated by 
looking at the particles’ steady state velocity distribution 
P v and at the corresponding static interparticle spacing 
distribution (at zero velocity) P r , both shown in Fig. 6 for 
increasing 0. The important points here are that small 
interparticle distances are energetically costly, and that a 
large spread of velocities relative to the center-of-mass de¬ 
notes larger frictional dissipation, according to the RHS 
of Eq. (3). At perfect alignment, short distances (col¬ 
ored column) are very frequent, which is energetically 
costly. At the same time the spread of velocities is mod¬ 
erate and so is friction. In the optimally misaligned case 
6> op t instead, the shortest distance becomes less frequent, 
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FIG. 6. (Color online) (a-c) The distribution P r of nearest- 
neighbor distances in the static relaxed configurations (Fd = 
0, Uo = 0.27) at rotation angles 0 — 0°, 0 = 6 opt = 4°, 
and 6 — 10°. (d-f) Corresponding velocity distribution P v of 
particles sliding under an external force Fd — 0.36. 


thus reducing energy as we already know. At the same 
time however P v develops longer tails at lower and higher 
particle velocities, both of which increase friction. At 
0 > # 0 pt finally the velocity spread drops and so does 
friction, the monolayer sliding less and less affected by 
corrugation. 

In conclusion, colloid monolayers in an incommensu¬ 
rate optical lattice develop, in full equilibrium and with 
realistic parameters, a small-angle structural misalign¬ 
ment, quite evident in moire patterns such as those 
of Fig. 3. Upon forced sliding, this misalignment can 
considerably increase the sliding friction, directly ex¬ 
tractable from the colloid drift velocity in experiment, 
using Eq. (3), over the hypothetical aligned geometry. 

The present results and understanding naturally ex¬ 
trapolate to the sliding of misaligned incommensurate 
lattices in contact such as, for example, physisorbed rare- 
gas or molecular submonolayer islands [6-8]. An inter¬ 
esting side aspect is in this case that misalignment trans¬ 
forms the orientation angle of a physisorbed island, gen¬ 
erally assumed to be fixed, into a continuous and possibly 
dynamical variable. The inertial sliding friction of such 
islands determines the inverse slip time in quartz crys¬ 
tal microbalance experiments [9], whose data must, at 
least in some cases, embed the frictional enhancement 
caused by misalignment when present. Even though the 
time needed to diffuse-rotate a ^ 100 nm-size island may 
be exceedingly large, the orientation angle distribution 
of islands will usually, under either stationary or sliding 
conditions, be continuous rather than delta-function like. 
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The general unavailability of 2D lattice orientation of in¬ 
commensurate rare gas islands (as opposed to full mono- 
layers, whose epitaxy generally differs) poses at present 
an obstacle to the investigation of these effects, which 
must nonetheless be considered as generically present and 
effective. The possible local misalignment of incommen¬ 
surate 3D crystals in contact and its potential role in 
sliding friction is an even less explored, but interesting 
issue which remains open for future consideration. 
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by the Swiss National Science Foundation through a SIN- 
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I. MODEL DETAILS 


We describe the colloidal particles as classical point¬ 
like objects interacting via a repulsive Yukawa potential 

Upp (r) = 2 exp (—r/A D ), (1) 


where r is the interparticle distance, Q is the coupling 
strength, and Ad is the Debye screening length. We re¬ 
strict particle motion to two dimensions, where colloids 
at rest form a triangular lattice of spacing a c in the (x, y ) 
plane. The periodic triangular corrugation potential 


W (r) 



U 0 w( r) 


27 TX 

2 cos-cos 

ai 


2?r y 

V3ai 


+ COS 


47 xy 

VS ai 


( 2 ) 


has strength Uo and periodicity a/, mimicking the exper¬ 
imental laser-induced optical lattice. The ratio p=ai/a c 
defines their relative commensurability. Here we focus 
on mismatched cases only, namely underdense incom¬ 
mensurate UI (p < 1) and overdense incommensurate 
OI (p 1). A detailed discussion of the approximations 
behind this model can be found in Ref. 2. 

Simulations are carried out using either constant-area 
periodic boundary conditions (PBC), or open boundary 
conditions (OBC) with a confining potential. In PBC, 
the purely repulsive particles are naturally confined and 
the total potential energy of N p particles is 


N v 




i= 1 


wi^ + l^Uppinj) 




(3) 


We use a triangular-lattice supercell defined by primi¬ 
tive vectors ai=L(l,0), a 2 =I/(0.5, VS/2) initially filled 
with a regularly spaced triangular monolayer of lattice 
constant a c = 1, taken as our length unit. The super¬ 
cell edge L and the corrugation potential wavelength are 
then chosen to achieve a given overall periodicity of the 
colloidal lattice and the periodic potential. After a num¬ 
ber of trials, we settled on simulations using L = 89 a c 
(N p = 7,891) and p = 89/96 ~ 0.927, corresponding to 
the fifth approximant in the continuous fraction expan¬ 
sion of 3/(1 + VS). 

In OBC a large square supercell of size L = 500 a c is 
partly filled by a circular island of N p ~ 30, 000 particles, 



FIG. 1. Average nearest neighbor spacing as a function of 
the distance from the center of the model colloidal island con¬ 
taining N p = 28, 837 particles, optimized at rest in OBC at 
T = 0, Uo = 0, and p ~ 0.927. The kinks at large radii are 
connected with the presence of line defects appearing in the 
relaxed configuration as a result of the growing inhomogene¬ 
ity. 


moving in the total external potential 

V x t(r) = G(r) [-Aa + W(r)], (4) 

where 

G(r) = exp ('ft)’ (5) 

is an unnormalized confining Gaussian of large width 
ac- As suggested by experiment [1] both the confin¬ 
ing amplitude Aq and the corrugation potential W(r) 
are controlled by the same Gaussian modulation. By 
fixing Q = 10 13 , A d = 0.03, A G = 1200, = 1200, 

the Gaussian confinement and the interparticle repul¬ 
sion have been balanced so as to yield an interparticle 
distance a c ~ 0.983 at the island center. As shown in 
Fig. 1, the lattice spacing is fairly constant in the central 
region and increases by a few percent toward the island 
periphery. The total potential energy has the same ex¬ 
pression as Eq. (3), but with V ext (ri ) in place of FF(r^). 
We carry out simulations at different lattice mismatch 
p = 0.927, 1.08, 0.84 set by choosing appropriate values 
of ai given the central colloidal crystal spacing a c . 

The equation of motion for the j-th particle displace- 
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Model Expression 

Typical Value 

Viscous coefficient 77 

6.3 x 10 -8 kg/s 

Length a c 

5.7 /im 

Force F = Qe~ a ‘ /X °/(Xoa c ) 

20 fN 

Energy E = Fa c 

1.1 x 10 -19 J 

Velocity V—F/r\ 

0.3 /im/s 

Power P—F 2 /rj 

6.3 x 10 - 21 W 

Time to = r]a c /F 

18 s 

Mass 771 — rfa c /F 

1.1 x 10 -6 kg 


TABLE I. Basic units for various quantities in our model and 
typical values mimicking the setup of Ref. 1. 


ment r j is 


mr j +r](r j -v d ) = -V rj 


y ] Cpp (vjj ) + Vxt ( ) 


(6) 


where is the drift velocity, giving rise to a Stokes’ drag 
force Fd=rjVd, experienced by each particle [2]. The over- 
damped dynamics typical of this system is achieved by 
adopting a sufficiently large value of the damping coeffi¬ 
cient rj. We use 77 = 28 in all simulations. All results are 
expressed in terms of the basic units defined in Table I, 
where values are inspired by those in the experiment by 
Bohlein [ 1 ]. Figure 3 in the paper refers however to a 
different experimental setup [4] corresponding to a softer 
colloidal crystal at a larger lattice spacing a c ~ 6.5 fi m. 
In that case a more appropriate value for the basic energy 
unit is E ~ &BT room = 4.1 x 10 -21 J. 


II. SIMULATION PROTOCOLS - STATICS 

In PBC the equilibrium configuration at each individ¬ 
ual value of U 0 is generated as follows. We start with 
a perfect colloidal crystal at 0 = 0 °, as sketched in 
Fig. la of the paper. A first series of zero-temperature 
(T = 0) damped molecular-dynamics simulations is car¬ 
ried out increasing the substrate potential amplitude in 
steps A U = 0.01 from Uo = 0 -A 0.5, relaxing the 
monolayer at each step until no appreciable movement 
of the colloids is observed. A second series of configu¬ 
rations is generated starting with a corrugation ampli¬ 
tude Uo = 0.5, decreasing it in steps A U = —0.05 down 
to Uo = 0 and performing an annealing protocol from 
T = 1 —»• 0 at each step. The stability of each of these 
annealed structures is subsequently tested by further car¬ 
rying out parallel T = 0 relaxations varying the corruga¬ 
tion amplitude in steps AU = ± 0.01 (both upward and 
downward), until Uo = 0 or Uo = 0.5 are reached. For 
each value of Uo, among all the configurations generated 
from these two procedures we eventually select the one 
with the lowest total energy. The resulting lowest-energy 



FIG. 2. A small region of the particle system optimized 
at rest in PBC. Particles are colored according to the local 
substrate potential W(x,y). (a) Aligned (9 — 0°) configu¬ 
ration obtained at Uo = 0.09. (b) Configuration obtained 

at Uo = 0.18, showing a spontaneous misalignment angle 
6 ~ 2.3°, whose effect is greatly magnified in the correspond¬ 
ing moire pattern, rotated by « 30° with respect to (a). An 
extended defect (at the left side) and a vacancy (pointed at 
by the arrow) appear, adapting the misaligned particle lattice 
to the constraint introduced by the PBC. 


configuration is usually aligned for small corrugation Uo, 
and develops spontaneous misalignment and defects for 
large enough Uo, see Fig. 2 . 

In OBC, equilibrium configurations are generated 
starting from a circular island cut out of a perfect tri¬ 
angular crystal (a c = 1 ) and fully relaxed inside the con¬ 
fining potential (Uo = 0 ). The density profile resulting 
from the Uo = 0 relaxation is drawn in Fig. 1. We then 
carry out damped-dynamics simulations at T = 0 in¬ 
creasing the corrugation Uo in steps AU = 0 . 02 . We 
follow the same protocol starting from islands rotated at 
different initial misalignment angles 0 ^, ranging between 
0° < 0i < 30°. From this procedure, we obtain energy 
curves like those in Fig. 2 of the paper, whose minima de¬ 
termine the optimal static configurations and equilibrium 
misalignment angle 0 op t- Note that in this process the fi¬ 
nal orientation 0 / of the island at the end of optimization 
is not identical to the initial angle 0 ^; however, it usually 
remains very close. Results are then given as a function 
of 0 = 9i. All relaxed configurations with 0 7 ^ 0 opt rep¬ 
resent locally stable (and generally metastable) states of 
the system. 


III. SIMULATION PROTOCOLS - SLIDING 
DYNAMICS 

In real sliding experiments the driving force is the 
viscous drag exerted on each colloid by a slow fixed- 
amplitude oscillatory motion of the experimental cell. 
The slow time scale yields long time intervals with con¬ 
stant drag force. Our dynamical simulations are there¬ 
fore carried out by applying to each particle the same 
external force Fj = r\Vd , usually along the x axis. F'd is 
kept constant during a finite time tp inversely propor¬ 
tional to the force value Fd itself (and thus the overall 
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FIG. 3. Static structure optimization. Total energy per 
particle as a function of 0 obtained in OBC at the overdense 
p = 1.08, for different corrugations Uq. The dashed line in¬ 
dicates the expected value 6>nm — 2.6° from Eq. (1) of the 
paper. 


sliding speed). After this time, the force sign is reversed 
for the same time duration. The particle motion being 
overdamped, there are no effects of inertia. The product 
Ffitp is chosen so that an isolated unconfined particle 
would move by a few lattice spacings Fdtp/v ~ 2 — 3 a c 
typically. 

When we study the sliding of the monolayer in OBC, 
the Gaussian confining potential is kept fixed in time. 
Only the homogeneous central part of the island is slid¬ 
ing, whereas the boundary tends to remain pinned dur¬ 
ing the slow oscillatory cycle. In order to exclude unde¬ 
sired edge effects, we select a square central region of size 
80 x 80 a^, containing ^ 7,500 particles. Aiming at de¬ 
scribing steady-state sliding, we discard an initial tran¬ 
sient of approximately 30% of the simulation time and 
compute averages considering particles which at t = 0 
were residing in this central block. 

Following Ref. 2 , the dissipated friction power is cal¬ 
culated according to Eq. (3) in the paper. Curves re¬ 
ported in Fig. 5 of the paper are obtained by applying 
this procedure to each of the misaligned static configura¬ 
tions produced following the protocol described earlier. 
We verified that due to the relatively large size of the 
island and the short duration of the simulations the ini¬ 
tial angular orientation of the island is preserved during 
the sliding trajectory (even at finite temperature). We 
carried out mainly T — 0 simulations integrating the 
equations of motions using an adaptive Runge-Kutta al¬ 
gorithm. When needed, thermal (Brownian) motion is 
simulated by adding suitable random Gaussian forces in 
a Langevin approach. 


IV. OVERDENSE INCOMMENSURATE 

Figure 3 displays the total energy as a function of 6 
for the overdense p = 1.08. The theoretical NM angle is 
#nm = 2.6°. The minimum is found in correspondence 
of 6 > op t — 3°, closer to the expected value than what was 
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FIG. 4. Sliding dynamics in OBC. Dissipated friction power 
Pfric as a function of 6 obtained at p — 1.08, external forces 
Fd = 0.18, 0.27, 0.36, and corrugation Uo = 0.18 (T = 0). 


found in the UI regimes. In this overdense case, moving 
toward the periphery the average colloid spacing tends to 
get commensurate to the corrugation potential, favour¬ 
ing smaller epitaxial angles. In the OI regime the misfit 
stress is confined in local compressions rather than di¬ 
lations. Due to the exponential form of the repulsive 
interaction the cost of compressions is higher, so that, 
at equal corrugation Uo, distortions are weaker than in 
the UI case. Consequently the dissipative sliding of the 
aligned and optimally tilted configuration is not very dif¬ 
ferent and the difference in pf r i c is reduced compared to 
the underdense regime. Nonetheless, far away from 6 opt 
a significant drop of the dynamic friction is still observed, 
as shown in Fig. 4. 


V. TEMPERATURE EFFECTS 

The NM theory [3] assumes zero temperature. With¬ 
out attempting a rigorous investigation of the problem 
here we anticipate, based on Eq. (1) of the paper and on 
finite-temperature simulations, that thermal effects do 
not affect the main conclusions obtained at T = 0. 

Within a mean-field theory, such as a quasi-harmonic 
approximation, one may still use T = 0 results with ther¬ 
mally renormalized constants. The two parameters defin¬ 
ing #nm are the mismatch ratio p = ai/a c and the sound 
velocity ratio 7 = cp/cp. The mismatch ratio is un¬ 
der control as one can always set the desired incommen- 
suration by fixing the average lattice spacing (a c ). Of 
course a c will be characterized by a broader distribution 
at T % 0. Changes in the sound velocities are also irrele¬ 
vant if the effects of temperature on the crystal properties 
act simply as a renormalization of the 2 -body interaction. 
Under this hypothesis 7 is expected to be weakly affected 
since, as it turns out, its value is mainly determined by 
the symmetries of the lattice [5]. Therefore in a mean- 
field treatment, the rotated epitaxy is expected to survive 
with only quantitative changes related to lattice thermal 
expansion. 

The peak in the dynamic friction is connected with 
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the increased interdigitation of the monolayer within the 
substrate occurring at 0 opt . At the optimal misalignment 
6> op t, solitons are enhanced, which is reflected in the larger 
distortions observed in the colloidal crystal (Figs. 4 and 
6b of the paper), and their motion is more dissipative. At 
room temperature (&BT room ~ 0.04 in simulation units) 
we observed that such solitonic superstructures are pre¬ 
served for corrugations Uo > 0.1 and dynamic friction 
does indeed increase near the optimal misalignement, as 
shown in Fig. 5. 



FIG. 5. The dissipated friction power pf r i C as a function 
of 6 at Troom- Simulations have been performed in OBC, 
fixing p — 0.927, external force Fd = 0.36 and corrugations 
Uo = 0.18, 0.27. 
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